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ABSTRACT 

Using the cosmological baryonic accretion rate and normal star formation (SF) efficiencies, we 
present a very simple model for star-forming galaxies (SFGs) that accounts for the mass and redshift 
dependences of the SFR-Mass and TuUy-Fisher (TF) relations from z ~ 2 to the present. The 
time evolution follows from the fact that each modeled galaxy approaches a steady state where the 
SFR follows the (net) cold gas accretion rate. The key feature of the model is a halo mass floor 



10 Mq below which accretion is quenched in order to simultaneously account for the 



observed slopes of the SFR-Mass and TF relations. The same successes cannot be achieved via 
a star-formation threshold (or delay) nor by varying the SF efficiency or the feedback efficiency. 
Combined with the mass ceiling for cold accretion due to virial shock heating, the mass floor Mmin 
explains galaxy "downsizing" , where more massive galaxies formed earlier and over a shorter period of 
time. It turns out that the model also accounts for the observed galactic baryon and gas fractions as 
a function of mass and time, and the cosmic SFR density, which are all resulting from the mass floor 
A^rnin- The model helps to understand that it is the cosmological decline of accretion rate that drives 
the decrease of cosmic SFR density between z ~ 2 and z = and the rise of the cosmic SFR density 
from z ~ 6 to z ~ 2 allows us to put a constraint on our main parameter Mmin — 10^^ Mq. Among 
the physical mechanisms that could be responsible for the mass floor, we view that photo- ionization 
feedback (from first in-situ hot stars) lowering the cooling efficiency is likely to play a large role. 
Subject headings: cosmology: observations — galaxies: high-redshift — galaxies: evolution 



1. INTRODUCTION 

To a good first order approximation, galaxies 
are either blue and active or red and passive, as 
indicated by the color bi-modality. In the past few 
years, it has been realized that blue star-forming 
galaxies (SFGs) lie on a tight relationship between 
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elevated SFRs are a natural outcome of merger-driven 
starbursts. However, the tightness of this relation, with 
rms scatter of less than 0.3 dex, indicates that SFR is 
not driven by merger-induced starbursts but rather by 
a continuous mass-dependent process that is gradually 
declining with time. 

The mean mass dependence and time evolution of 
the relation between SFR and stellar mass can be 
summarized by the expression 



Af* = 150Af/ii (1 



^3.2 



yr" 



(1) 



SFR-Afj, relationship is analogous to the red sequence 
for passively evolving galaxies and is sometimes referred 
to as the SFR sequence. Every multi-wavelength survey 
has shown that the specific SFR (sSF R=SFR/M^) is 
higher for l ower mass SFGs (e.g [Brinchm ann et al.ll2004l: 
Bell et al.! 120051: INoeske et al.l [2007at iDrorv fc AlvarezI 
2008; Oli ver et all I2010D . i.e.. the mass index of the 
SFR-Af* relationship is less than unity. The SFR 
sequence has evolved by a factor 20 at a given stellar 
mass from z ^ 2 to the present time. This strong 
evolution of the SFR sequence implies that distant 
z - 2 SFGs wit h > IQ'^°-^Mq h a d SFRs in excess 
of 100 Mp, yr -^(Shaplcv et al.' 2003; Er b et al.l [2006al: 
IGrazian et all 2007: Daddi et al. 2007). Locally, such 
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where Af,,ii = M^/W^Mq, (1 + z)3.2 = (1 + z)/i.2, 
p ~ 0.8, and q ~ 2.7 in the redshift range z — 
— 2.0 Reproducing the characteristic mass and time 
dependencies of the SFR seq uence is a challenge for 
models of galaxy formation (e.g. lDavel2008l : lDamen et all 
12009a'). 

Equation [1] is important as it is very reminiscent of 
the halo mean growth rate, which has been shown to 
be Afh (X Afh'' (l + zY with a mass in dex s greater than 
unity s - 1.1 |Neistein fc De kcl 2008. iGenel erall[200l 
IMcBride et al.112 009') that is set by t he shape of the initial 
dark matter (DM ) power spectrum ( Neistein et al.|[2"006l : 
iBirnboim et al.l [20071 ). The similarity motivates a closer 
investigation, and we will show that, indeed, there is a 
strong intimate connection between the growth of halos 
and the SFR sequence. 

Another relevant scaling relatio n for (disky) SFGs is 
the well-known TuUy-Fisher (e.g. iTullv fc Fished 1 19771) 

® There arc marginal indications for a variation of p from near 
0.7 at 2 = - 1 (Brinchmann ct al. 2004; Noeske ct al. 20073) 



to about 0.9 at z '~- 
IPannella eran[2009l) . 
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relation, which correlates stellar mass Mi, and maximum 
circular velocity Vi„ax- The TF relation a ppears to be 
already in place at high-redshifts z > 1 (iKassin et al " 



20071 iPuech et al.]l200a lEpinat et al1l2009HCresci et al 



20091) and has evolved by a factor of 2.5 at a given mass 
(jCresci et al.,i2009.) from z ~ 2.2 to the present. It is of 



the form: 



(2) 



where m 4 to a 10 % accuracy locally (e.g. 

ITuUv fc Piercg[2000l:lMcGau gh 2005; Meyer et al. 2008). 
As for the SFR sequence, Equation [5] is reminiscent of 
the virial relation for DM halos, namely, Mh oc V^. 
Thus, both the SFR sequence and the TF relation behave 
similarly to their DM counterparts, but have slightly 
different mass indices, i.e., they are tilted with respect 
to their DM counterparts. 

Aside from our knowledge on these global properties, 
our detailed understanding of individual SFGs has 
also improved greatly thanks to spatially resolved 
kinematic studies of the ionized gas in z ~ 2 SFGs 
'e.g. SINS survey. iFor ster Schreiber et al.' 

Hoo 



Genzel et al I [2006. 
Cresci et al.M2009D . 



'2006', '2009; 



12008; Shapiro ct al. 2008, 2009; 
The SINS survey, consisting of 80 
revealed that 30%-50% of SFGs 
"^^■^ are gas-rich thick rotating 



z = 2 SFGs, has 

with ~ 10^°-^" 

disks (see also Ivan Starkenburg et al.ll2008l : IWright et al.l 
I2007I) . with the rest either dominated by dispersion 
veloci ties or showing t h e sign a ture of r ecent mergers 
(e.g. iLaw et al.l 120071 120091: IForster Schreiber et al.l 
l2009l) . The high-redshift SFG disks are different 
from local spirals in many ways. In particular, 
they are t hick and correspondingly of high velocity 
dispersion ("Forster Sc hreiber et al.l 120061 : IGenzel et al.l 
[2008; Cresci ct al. 20 0j), and they are often made of 
giant SF 'clumps ' ICowie et al.| [l995l: IGenzel et all[200l 
lElmegreen et al.l 120071 Forster Schreiber, Shapley et 
al. in prep.). This feature indicates wild gravitational 
instability, which is a likely outcome of a high gas density 
jNoguchi 1998; Immch ct al. 2004; Bournaud et al. 2007; 
lEhnegreen et al.l 120081: IDekel et al.l 12009b V Indeed, 
direct evidence for high gas fractions at high redshifts 
is now mounting jEvh et al.l l2006bl: iDaddi et al.l 120081: 



iTacconi et alll20ld iDaddi et al.ll2010al) thanks to rapid 
progress in the sensitivity of mm interferometers. 

The global scaling relations and the recent surveys of 
galaxy kinematics raise several outstanding questions: 
(1) What causes the high SFRs at z ~ 2? (2) What drives 
the evolution of the SFR sequence and TF relation from 
z 2 to the present? (3) What drives the cosmological 
evolution of the average SFR density? (4) Why did 
more massive galaxies form their stars before less massive 
SFGs? (5) Why are z = 2 SFGs so gas rich? 

In this paper, we address these questions in the context 
of the cosmological growth of DM halos. In particular, 
we construct a very simple model that ties the SFR 
sequence and TF relation to their DM counterparts. The 
important feature of the model is the suppression of 
the accretion below a mass floor at M^i^ ~ 10^^ Mq, 
which, as we will show, has several other important 
consequences. We focus our attention at z ~ 2, i.e., 
when the universe was just 3 Gyr old. 

This paper is organized as follows. In Section [2l we 
present the model. In Section |3l we show that the SFGs 



reach a quasi-steady state and demonstrate the impact 
of the mass floor Mmin on the SFR sequence and the 
TF relation. In Section |4l we find that the mass floor 
naturally delays the star-formation activity and leads to 
downsizing. The model also simultaneously accounts for 
the baryonic and gas fractions as a function of mass and 
time. In Section [5l we put direct constraints on the 
numerical value of Mmin from the cosmological history 
of star-formation density. In Section |6l we discuss the 
model limitations and present our conclusions. Finally, 
in Section [71 we discuss the possible origin of the mass 
floor. We use throughout the standard ACDM cosmology 
with the parameters f2,n = 0.3, Ha = 0.7, h = 0.7 and 
0-8 = 0.8. 

2. THE RESERVOIR MODEL 

In this Section we present our 'reservoir' model and 
its two major ingredients (Section 12. ip . namely, the 
accretion efficiency, and the SF efficiency. We discuss the 
mass and redshift dependences of these two parameters in 
Section [221 The model ingredients are then summarized 
in Section [2.31 We will use this model as a learning tool 
to gain insights on the role played by several key physical 
parameters. 

2.1. The basic equation 

We consider a galaxy with its dark halo as a reservoir 
that is fed by a source and is emptied into a drain. The 
source represents the amount of newly accreted cold gas, 
and the drain is the gas consumption into stars as well 
as outfiows. The basic equation of our model is the 
differential equation expressing the conservation of gas 
mass, 



-AL 



(1 - R)M, - M, 



gas, out ; 



(3) 



where Mgas,in is the gas accretion rate, (1 — R)Mi, is the 
net SFR corrected for the recycled fraction, and Mgas.out 
is the mass outflow rate. For our purposes, i? is a kept 
time-independent given its slow de pendence on stellar 
popu lation age T for T > 10^ Gyr (jBruzual fc CharlotI 
120031) ■ i.e., at z < 4. We discuss the limitations of this 
assumption in Appendix A. 

To be as general as possible, we included an additional 
drain in our model (Equation [3]) representing any 
outflowing gas. Observationally, the outflow rate 
-^^gas.out is observed to be roughly proportional to the 
SFR Mgas.out = a X SFR, with a c± 1.0 (iHeckman et al.l 
[2000t IMartinI [20051 IRupke et al.l[200l lErbl 120081) . "V\^ 
note, however, that some fraction of the supernova (SN) 
driven winds is likely to be recycled, especially in the 
most massive halos, where much of the outflowing gas 
falls back into the galaxy and can boost the stellar mass 
at z = by 33% (simulations bv lOppenheimer fc Davel 

I2008D . Hence, the net outflow rate Mgas,out may be less 
than the SFR, i.e., Mgas,out ^ SFR. 

Regardless of the numerical value of the 
proportionality constant a between outflow rate 
and SFR, Equation [3l can be re-written as 



(4) 



where a includes the corrections for recycling and 
outflows. This equation expresses the trivial fact that 
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the gas reservoir Afgas will be filled up or get emptied 
depending on the relative power of the source and the 
drain terms. As discussed in section [3. 11 the system will 
self-regulate itself to a steady state (Afgas — 0) where 
SFR is proportional to the cosmological accretion rate. 

We stress that the feedback term will have little impact 
on the scaling relations, such as the SFR sequence, or 
the sSFR. Indeed, if half of the accreting gas is entrained 
and expelled by the SN-driven winds, the SFR will be 
lowered by a factor of 2. As a result, if the SFR is lower 
by a factor of 2, the stellar mass Af^ = / dt SFR(t) will 
be lower by the same factor, and the SFR sequence will 
remain unchanged. 

2.1.1. Halo growth and gas accretion 

The gas replenishment term Mgas,in in Equation |3] 
is require d observationally. F o r example, the G-dwarf 
problem (Ivan den Berghl Il96l ISchmidtl MM) c alls for 
a sig nificant amount of newly accreted gas (e.g. iLarsonl 
Il974| ). Furthermore, the gas depleti on timescale in local 
massive galaxies is ~ few Gyr (e.g. iWong fc Blitz! [20021 : 
iJames et al.l I2008D . i.e., much shorter than the time 
required to build their stellar masses. Similarly, for z ^ 2 
SFGs, the gas consumption timescale of less than 0.5 Gyr 
is shorter than their typical s tellar ages of 1-2 Gyr (e.g. 
lE7bl[200llTacconi et al.llMol ). thus requiring intense gas 
accretion. 

Cold gas accretion is also required theoretically as 
it is a natural consequence of the 'cold-acc retion' 
regime (jBirnboim fc Dckel"2003': Keres et al.ll2005f ) when 
the cooling time is s horter than the dynamical time 
(jWhite fc Frenld 119911 ). Furthermore, since the cold 
gas is not shock-heated to the virial temperature, high 
gas accretion efficiency is a natural outcome of efficient 
penetration of cosmological cold streams into the inner 
galaxies at high redshif t, as seen in hydrodynamical 
cosmo logical simulations (jOcvirk et al.ll2008t iDekel et"all 
l20n9aD . 

In this context, it is the halo growth rate of DM 
halos that is regulating the baryonic accretion. The 
halo growth rate is by now well understood based on 
A^-body simulations and the exten ded Press-Sc h echter 

(EPS) analyt ic fo rmalism (lEfsta thiou e tlUI 119851: 

Wechsler et "alll2002t Fran den Boschl l2002: Spr ingel et al] 



2006[) . An approximation for the average mass growth 



rate of DM halos of virial mass Afh at redshift z in 
cosmologica l A^-body simulations (van d en Bosch|[2002l: 
iGenel et al.l [2008: McBride et al. 2009), which is also 
understood using t he extended Press-Sch echter (EPS) 
analytic formalism ([Neistein fc Deke'ill2008| ). is obtained 
by the fitting function 



Afh ~ 510Mh^i2(l 



2^)3.2 yr 



(5) 



where Mh,i2 = Mh/lO^^ Mq, (1 + 2)3.2 = (1 + ^)/3.2, 
t ~ 2.2, and s ~ 1 .1, with the estin iates f or s ranging 
from 1.08 to 1.14 ([Neistein fc Deke]| [200l IGenel et al.l 
[20081 IMcBride et al.ll2009l) 

Given the average halo growth rate, the corresponding 
average gas accretion rate is 

Mgas,in = fin /b Afh 

^ 90 ein /b.0.18 A4;i2 (1 + Mq yr-1 , (6) 



where /b,i8 = /b/0.18 is the cosmic baryonic fraction, 
and ^ 1 is the accretion efficiency expressing the 
effective fraction of the baryonic matter that is actually 
accreted as cold gas into the galaxy. 

The accretion efficiency ei„ is a very important 
parameter and must be within the range 0.5-1.0 
for several reasons. Theoreti cally, is hig h 



under the cold- accretion regime (jDekel et al.l 12009a ). 
Observationally, this efficiency is supported by the fact 
that the observed SFR in massive galaxies at z ^ 2 
matches within a factor of tw o the maximum predicte d 
gas accretion rate (IGenel et a l. 2008; Dekel et al.ll2009al) . 
In addition, IGenel et al.l ( 2008i ) showed that the halo 
major merger fraction in the mass range of interest 
logAfh = 12 is low and consistent wit h the low 
merg er fractions in the SINS survey (S hapiro et al.l 
I2008t lF5rster Schreiber et al1l2009[ ). provided that the 
accretion efficiency ein is grea ter than 50% (see F i g. 1 o f 
IGenel et al.ll2008D . Recently. [Bauermeister et al.l (|2010f ) 
compared the inferred 'external' gas accretion to the 
DM accretion rate and found this efficiency needs to be 
> 70%. In light of all these results, we will adopt a 
fiducial value of ein = 0.7. 

As mentioned in the Introduction, this baryonic growth 
rate (Equation [6]) is very reminiscent of the SFR-Mass 
relationship (Equation [Ij , given the similar near-linear 
increase of the accretion rate with mass. However, the 
mass index s in the accretion rate is slightly larger than 
unity, while the mass index p in the SFR sequence is 
somewhat smaller then unity. The clear implication 
from this difference is that the specific accretion rate 
increases with mass for the DM component, while the 
sSFR decreases with mass. Thus, if there is a connection 
between Equation [6| and Equation [l] the difference has 
to be explained (see Section [3]). 

2.1.2. Star formation rate 

Star formation, a very complex, local and inefficient 
process (giant molecular clouds (GMCs) typically 
turning 1-2% of the gas mass into stars in a free-fall time) 
is the primary drain of the reservoir (in Equation [3]). 
On galaxy scales, the amount of gas consumed is well 
described by the empirical relation: 



5^SFR — fsfrSgas/idyn, 1-6. 
SFR — Esfr Afgas/ (^dyn) , 



(7) 

where idyn = fii/2/K; is the galaxy dynamical time, 
and Csfr is the SFR efficiency parameter. For a 
marginally unstable disk with T oomre Q p arameter 
Q 1, it can be shown ( Martin fc K ennicutH 

120011: iKrumholz fc Thompsod |2007|) that Eq^iation [7| is 
equivale nt to the trad. i tional Kennicutt - Schmidt (KS) 
relation (jSchmidti [l95l iKennicuttI [19981) . which relates 
the surface densities of gas and SFR via Ssfr ^g as- 

The locally inferred value by lKennicut3 (|1998[ ) for 
the SFR efficiency is Csfr — 0.02 and the most recent 
advancements in constraining the KS rel ation at z ~ 
2 sh o w no eviden c e for any evolution (jGenzel et al.l 
[20Tot iDaddi et all l2010b[ ). In addition, at z ~ 
2, we know that massive disky SFGs, with M^, > 
1O^°-^M0, extend to half-light radii of R1/2 ~ 4 kpc 
(iBou che et al.l l2007al ) and the y rotate with circular 
velocities Vc ~200 km s~^ (jForster Schreiber et al.l 
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120061) . Expressing the quantities in these units, the 
orbital time can be written as 

At z > 2, we assume that the orbital time scales with the 
halo dynamical time, Rvir/Vh, which is a fixed fraction, 
about 18%, of the Hubble time at the given redshift, thus 
proportional to (1 + z)~^^^. 

We note that our results on the scaling relations 
presented in 13.21 are completely independent of the SF 
efficiency (i.e., egfr and idyn)- This is because the SFR is 
going to be driven by the cosmological accretion rate, as 
explained in Section IXTl 

2.2. Mass dependence 

The two efficiency parameters egfr and ein appearing 
in Equations |6]-[7] may in principle vary with mass and 
redshift. The SF efficiency Csfr does not appear to evolve 
with redshift. Indeed, the KS relation seems to be in 
place as early as z = 2 and shows no significant evolution 
between z = and z = 2 (jBouche et al.ll2007bD . 

On the other hand, the accretion efficienc y e\„ must 
be a strong function of halo mass (e.g. White fc FrenkI 
I199H ). At the massive end Mh ^ M^ax ^ IO^^Mq, and 
at low redshifts, the cold flows are not expected to be in 
a form of narrow, dense streams and they therefore fail 
to penetrate through the shock-heated halo gas. This 
may play a significant role in quenching star formation 
and causing t he transition of b l ue galaxies onto the 
red s equence (jPekel fc Birnboiml 120061 : iCattaneo et al.l 
|2006( ). To model this effect, we set a ceiling for cold-gas 
accretion at all redshifts Q 

ei„(Afh) = if A4 > Mmax. (9) 

At the low-mass end, there is strong evid ence for a drop 
in th e efficiency of galaxy formation (e.e. iShankar et alJ 
2006'; 'van den Bosch e t all 120071: iBaldrv e t all 120081: 
Kravtsov 2010; Mostcr et al.ll2010l: |Guo et al.||2010D. In 

the spirit of our simple toy model approach, we make the 
ansatz, 

ei„(Mh) =0 if Afh < M„i„, (10) 

i.e., we set a sharp mass floor for ein(Mh) below Mj^in- 
For any reasonable physical mechanism that may be 
responsible for the suppression of accretion and SFR in 
low- mass halos (see Section [71), modeling the effect as a 
sharp cutoff is clearly a very crude approximation, but it 
is useful as a first attempt in capturing the key features 
with a minimum number of parameters. 

This accretion floor means that the effective accretion 
of cold gas is suppressed in halos less massive than 
-^min- Halos morc massive than A/min accrete baryons 
that are embedded in other merging halos, which could 
themselves be more massive or less massive than Afmin, as 
well a s smooth gas unbo und to any halo (jStewart et al.l 
120081: iGenel et al.l I2O1O0 . The relative contribution of 
these different contributions leaves sufflcient margins to 
meet an accretion efficiency of 0.7. 

''' Cold accretion can occur above this mass threshold at redshifts 
beyond z = 2 when the filament cross-section is much smaller than 
the virial radius. However, the number density of such halos with 
mass 10^^ (and above) is low. 



The accretion efficiency ein at a given mass is likely to 
decline in time, as the overall fractions of stars and hot 
gas {'S> lO^K) in the inter-galactic medium (IGM) grow. 
We model this by incorporating a redshift dependence 
ein(z) = Cin f{z), to allow for a decrease in cold accretion 
efficiency from z = 2 to z = 0. For simplicity reasons, 
we use a function linear in time with the boundary 
conditions f{z = 2.2) = 1 and /(z = 0) = 0.5. Note 
that our main results on the slope of the scaling relations 
are independent of this assumption since a change in 
Ein changes both SFR and M^, by the same amount. 
However, as we discuss in Section 3.4, this affects the 
evolution of the SFR sequence. 

Figure [Tfleft) shows the mass growth histories for our 
modeled halos and the model accretion efficiency ein(A'/h) 
with the upper and lower cutoffs at Afmax and Afmin- 
Figure [Ij right) shows the maximum baryonic accretion 
rate fs Afh available for each halo as a function of 
time. More massive halos reach the Afmin threshold first 
and cross Afmax over a shorter time scale, leading to 
downsizing (see section ITT]) . 

2.3. Overview of model ingredients 

At this stage, it is worth summarizing our model and 
its parameters. The model parameters are as follows: 

1. The key assumption in our model is the halo mass 
floor for cold gas accretion Afmin, on the order of 
10^" - 10" Mq. As we wiU show in Sections [X^ 
andO Afmin ~10^^ Mq wfll be required to match 
observations. 

2. The halo mass ceiling for cold gas accretion 
-Mmax, set at a value of 1.5 x 10^"^ Mq to 
match th e characteristic thr e shold for virial shock 
heating (jPekel fc Birnboiml 120061 ) and for best 
reproductio n of the observed feat ures of the galaxy 
bimodality (jCattaneo et al.ll2006[ ). 

3. The cold gas accretion efficiency ein{Mh,z). As 
illustrated in Figure [U it vanishes below A^,„in and 
above Afmax, and is set to 0.7 at z > 2. Below 
z < 2, we use ein(z) = /(z) x 0.7, where /(z) is a 
function that is decreasing linearly with time with 
the boundary conditions /(2.2) = 1 and /(O) = 0.5. 

4. The SFR efficiency Csfr = 0.02 set from the 
observed KS relation, and the orbital timescale 
tdy„ = 2xl0^ yr(l + z)^i-5. 

Our fiducial model described above is the "accretion 
floor" model, dubbed hereafter 'accFloor', where the 
accretion of cold gas, and therefore star formation, are 
totally suppressed below A/min- We emphasize that our 
terminology 'accretion floor' means that the effective 
accretion of cold baryons is suppressed in halos with 
Afh(z) < Afniin- Either these baryons are prevented from 
cooling or they are prevented from entering and reaching 
the central object. 

In Section [3l we will analyze two radically different 
models for comparison (listed in Table [1]): (1) in the 
other extreme 'noMmin' model, the mass floor is not 
applied at all, i.e., we allow both gas accretion and star 
formation to proceed unperturbed below Afmin and (2) 
in the intermediate 'sfrFloor' model, we do allow gas 
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z z 

10 4 2.2 1.2 0.5 10 6 4 2.2 1.2 0.5 




IQ9 IQIO IQ9 IQIO 

Time [yr] Time [yr] 



Fig. 1. — Left: The growth history M^^{z) for each modeled halo as a function of redshift calculated according to Equation |5] Above 
the quenching mass Mmax, accretion is increasingly hot. Below the quenching mass Mmax, cold accretion brings in gas efficiently with 6-,^. 
Importantly, below a minimum mass Mmin, accretion is inefficient. Right: The maximum baryonic accretion rate inferred from the growth 
rate for each modeled halo. The final halo mass M^Jz = 0) is labeled. More massive halos reach the Mmin threshold first and cross A/max 
over a shorter time scale (downsizing; see Section l4.1l l. 



Model 


Afmax 


Afmin 


Floor 


Feedback 




Mq 


Mq 




a 


accFloor 


1.5 X 10^2 


10" 


accr. 





noMmin 


1.5 X 10^2 


n.a. 


n.a. 





sfrFloor 


1.5 X 10^2 


10" 


sfr. 





accFloor+ 


1.5 X 10^2 


10" 


accr. 


0.6 



TABLE 1 
Definition of our fiducial models. 



accretion below Mmin, but forbid star formation there. In 
this case, the cold gas accumulates until the halo becomes 
more massive than M„iin, and it then turns into stars. We 
also investigate the impact of the outflow term Mout in 
Equation [3] with a model we term 'accFloor+'. However, 
as described in Section[5?TJ this term will have no impact 
on the scaling relations. 

Figure [2] highlights the main differences between the 
three alternative models, 'noMmin' (left), 'accFloor' 
(middle), and 'sfrFloor' (right). In the top panels, 
we show the assumed variation of ein with redshift 
(thick solid lines, right axis) and the resultant evolution 
of the accretion rate (solid line), baryonic accretion 
rate (dashed line), and the SFR (dotted line) as a 
function of redshift. In the bottom panels, we show the 
resultant evolution of the DM halo mass (solid line), 
the maximum baryonic mass (dashed line), and the 
cold gas mass (dotted line). These three models differ 
significantly in their initial behavior. Compared to the 
'noMmin' model, the 'accFloor' and 'sfrFloor' models 
add a (mass-dependent) delay in the evolution of SFR. 
In contrast to the 'accFloor' model, the 'sfrFloor' model 
accumulates cold gas till the halo reaches Afmin- Both of 
these models lead to no SF below Mmin, but have rather 
distinct physical interpretations and consequences. 

2.4. Numerical approach 



Using this simple reservoir set up (Equations |3] and [7]), 
we solve for its gas and stellar content at each redshift as 
follows. We first construct average mass growth histories 
Mh(Mo, z) for the main progenitors of halos of masses in 
the rang e Mn = IQiO-^-i^ Mr, today, as in Figure [U 
following iNeistein fc Dek cT (2008") (see their Equations). 
We then numerically integrate Equations [SHZl for Afgas 
(and M-):) from z = 10 to any z until z = 0, with the 
initial conditions A/gas = Mi, = at z = 10. The 
stellar mass Af^, is updated at each time step according to 
NU = (1 - i?) / SFR dt, where R = 0.52 is the recycled 
gas fraction for a Chabrier initial mass function (IMF) 
(IBruzual fc Charlotl[200l Fl. 

3. RESULTS 

3.1. Steady state 

Before showing the numerical results, it is worth 
understanding the key behavior of the reservoir model 
from a qualitative inspection of Equations [3] and [71 
When Afgas in the galaxy is still low, the SFR is low 
by the KS law (Equation [7]) , and it is smaller than the 
accretion rate dictated by the cosmological environment. 
The gas reservoir then gradually fills up until the SFR 
becomes comparable to the accretion rate (see Figure [S]). 

* For our purposes, ij is a constant given its slow dependence on 
stellar population age T for age T > 10^ Gyr. 
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Fig. 2. — For a logA/h = 11.6 halo at z = 2.2, we show the behavior of each of our three models, namely when M^^^ = ('noMmin', 
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the same stellar mass Af* as the 'noMmin' model, given that the SF is simply delayed, and the amount of gas accreted remains the same. 
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Fig. 3. — For our fiducial model ('accFloor'), the bottom panel 
shows the SFR (solid line) and the maximum accretion rate 
(^in /b X Afh) (dashed line). The top panel shows that below 
z ~ 5.5, the accretion time scale (tacc = Afij/Afh) is longer than 
the SF time scale (igf^ = Afgas/SFR). The reservoir model reaches 
rapidly a quasi-steady state in which the SFR scales with the 
accretion rate. 



The galaxy enters a gwasi-steady state where the SFR 
is essentially set by the accretion rate, SFR~ Mgas,in- 
Would the gas reservoir be temporarily overfilled (for 
its current SFR), the SFR would then be larger than 
the accretion rate, and the galaxy will return to the 
quasi-steady state. 

Figure 13] demonstrates the quasi-steady state behavior. 
For our fiducial model ('accFloor'), the bottom panel 
shows the SFR (solid line) and the maximum accretion 
rate (/s Afh) (dashed line). The reservoir model reaches 
rapidly a gwasi-steady state in which the SFR scales with 
the accretion rate. 

Because the accretion rate varies with redshift, the 
steady state can only be achieved when the timescale 
associated with SFR is comparable to or shorter than 



the accretion timescale. The SFR timescale, using 
Equation [HI is 



gas 



(idy„,7/2) £3-^,0.02 (1 + zY^^" Gyr , (11) 

where torb,7 = idyn/lO"^ yr, and esfr,o.o2 = esfr/0.02. From 
Equation [SI the accretion timescale is 



tacc = Mh/A4 - 2 M^-°-' (1 + z)3:^ 'Gyr. 



(12) 



Therefore, the condition that the SFR timescale be 
shorter than the accretion timescale 



^sfr — 



(13) 



is met at redshifts z ^ 7. 

The steady-state solution SFRoc Mgas,in has several 
important implications. First, it implies that the SFR(z) 
is driven by the (net) accretion rate as illustrated 
on Figure [3] As a result, at around z ^ 2, SFR 
maintains a high and slowly varying value for a few 
Gyrs (Figure [TJb)), which means that the average SFR 
is expected to be comparable to the instantaneous SFR. 
This is indeed supported by the SINS sample at z ^ 2, 
where the birthrate parameter, defined as the ratio of 
instantaneou s to past SFR, i s estimated to be 6 ~ 1.2 
([Forster Schreiber et al.ll2OO90 . 

Second, the resulting SFR, gas and stellar masses are 
thus independent of the initial conditions given the rapid 
growth of the gas mass (see Figure [2]). 

Lastly, when the steady state condition (Equatior fTSl) 
is not met, the gas reservoir is being filled up at a rate 
faster than its consumption rate. In this case, SFR is 
time dependent, of the form i", and sSFR is proportional 
to 1/t. 

3.2. Effect o/Mmin on scaling relations 

Our goal is to better understand the global scaling 
relations of z ~ 2.2 SFGs, and the arguments presented 
in this Section are, however, applicable at z = 0, z = 1 
as well. Thus, we take the naive point of view that 
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Fig. 4. — The z = 2 SFR sequence predicted from our three fiducial models with no mass floor Af^jin (left), an accretion floor (middle), 
and an SF floor (right). The modeled points are shown as filled (open) squares, for halos below (above) the virial shock mass A/max- 
The observed SFR, sequence from PanncUa et al. (2003) and IDaddi et al. (2007) is shown as gray shaded areas. The dotted line in each 
panel shows a line of slope 1.1 expected from the global accretion rate (Equation [Sj . The 'accFloor' model provides the best match to the 
observed SFR sequence with SFRoc AfJ **. 
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V„„ [km/s] V___,^ [km/s] V_^,^ [km/s] 

Fig. 5. — The z = 2 TuUy-Fisher relation predicted from our simple model with no mass floor A/min (l^ft), with an ac cretion floor 
(middle), and with an SF floor (right). The calculation for a model with no A/jnin is shown as dashed lines. The data fromlCresci et all 
(|2009i ^ are shown as gray circles. The 'accFloor' model provides the best match to the observed TF relation with Af* oc V^^^ (solid line). 



the slopes of the scaling relations do not evolve with 
redshift and we stress that our aim is to demonstrate 
the key role played by the mass floor Afmin in connecting 
the baryonic scaling relations to their DM counterparts. 
Any evolution in the slopes of the SFR sequence or TF 
relation will lead to an evolution of this mass floor. 

Figure |4] shows the SFR sequence predicted by 
the three fiducial models ("noMmin", "accFloor" and 
"sfrFloor"), in comparison with the observed SFR 
sequence at that redshift from lPannella et al.l (|2009( ) and 
IDaddi et all ([2007). The data show clearly that the 
sSFR is higher for lower mass galaxies (the mass index 
p < 1), whereas the expectation from the DM relation 
would predict an inverted trend since its mass index 
s > 1. The slope s =1.1 associated with the predicted 
accretion rate as a function of halo mass is shown as the 
dashed line. Figure [S] compares the corresp onding z = 2 
TF re lation with the observed relation from iCresci et alj 
()2009() . Similarly, the data show a tilt from the slope 
3 associated with the virial relation Mh oc (dashed 
line) . 

These two figures highlight our first significant results: 
the 'accFloor' model provides the best match to the two 
scaling relations. The mass fioor A/min affects the slopes 



of both the SFR sequence and the TF relation in the 
right direction. It tilts the SFR sequence away from the 
1.1 dark accretion rate, and steepens the TF relation 
away from the 3 virial slope. This effect can easily be 
explained. 

When a mass floor is not applied, once the system 
reaches the steady state solution, SFR follows the halo 
accretion oc iWh and thus scales with halo mass 
according to Mi, oc M\^^^ . Given that the halo mass 
is Mh — /p A/h dt, and that the stellar mass is A/j,(t) oc 

Jp Mi, dt, stellar mass and halo mass are proportional to 
each other: M^, oc A/h- Hence, we expect in this case 
Mi, oc Af^-^, as seen in the left panel of Figure SI 

When a mass floor for accretion is imposed (model 
accFloor), accretion (hence SF) is suppressed until the 
time imin when the halo reaches Afmin, which is a 
decreasing function of ATh, and thus of Mi,. The stellar 
mass is now an integral from tmin(A'-/h) to t: Mi,{t) oc 
Mi, dt. The less massive systems had less time 
to grow stars and Mi, is now smaller at a given halo 
mass, by a larger factor for a smaller halo. This effect is 
responsible for the required tilt both in the SFR sequence 
and the TF relation as seen in the middle panels of Figs. 2] 
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and [51 

In the sfrFloor model, all the gas that has accumulated 
in the galaxy until t^in turns into stars after tmim so the 
total stellar mass after that time is expected to be the 
same as in model 'noMmin', as seen in the right panels 
of Figs, m and [SI 

In summary, a mass-dependent delay in the gas 
accretion is found to be necessary to better reproduce the 
observed slopes of the scaling relations. In order for the 
time delay to have a significant effect, we find that Mmin 
needs to be within 1 dex of Mi„ax, i-e., Mmin ~ 10^^ Mq. 
In Section [SJ we will put a stronger constraint on the 
numerical value of M^i^. The physical origin of the mass 
floor Mmin is addressed in Section [7] 

This mass-dependent delay can only be achieved in our 
model by the mass floor M^i^. Indeed, these results are 
completely independent of the SF efficiency and of 
the feedback term in Equation [3l A change in these two 
parameters will lead to a change in both and SFR, 
i.e., leaving the SFR sequence intact. A change in the 
accretion efRciency ein via its mass dependence is the 
only mechanism that we found to affect the slope of the 
scaling relations. 

3.3. Connections between the two relations 

It is very interesting to realize that the TF relation and 
SFR sequence are in fact tied to each other. Indeed, the 
steady-state solution oc Mh*, together with the SFR 
sequence Mi, oc implies the TF relation (A/^ oc V[™). 
Indeed, using the virial relation Afh oc V^, we find 

m — Zs/p . (14) 

With the indices being s ~ 1.1 and p ~ 0.8, the TF 
slope is indeed m 4 as observed (e.g. iMever et al.l 
I2008t lMcGaughll2005l ) . Note that this argument holds for 
various values of the indices s and p as long as s/p ~ 4/3. 
As noted in Section[Tl there are margi nal indications that 
jp varies from near 0.7atz = 0— 1 ([Brinchmann et all 
[2004: Nocskc ct al. 2 007al) to about 9 at z ~ 2 
TOaddi et al.l 120071 : iSantini et al.l 120091 iPannella etaP 
120091 ^. These subtleties do not change our arguments 
given that all observational evidence show that p < 1 
(sSFR higher at lower mass), while the DM counterpart 
has a mass index s > 1 (specific accretion rate higher 
at higher mass). Would this redshift dependence of the 
mass index p be confirmed, it would point toward a 
redshift evolution of Mmin . 

Given that the slopes of these two scaling relations are 
coupled together and that, as we will show, the accretion 
rate regulates the evolution of the SFR sequence, we 
expect the accretion rate to regulate the evolution of 
the TF relation as well. However, a more detailed 
analysis of the TF zero-point normalization requires 
knowledge of the explicit relation between the virial 
velocity of the halo, which we crudely referred to so far as 
Vhj and the observed maximum rotation velocity Vmax- 
Early models of galaxy formation assumed halo adiabatic 
contrac tion due to the infall of baryons into the halo 
center ([Blumenthal et al.l 119861: IMo et al.l [1991 . which 
gives typically Vmax — 1.6-1.8 Vh for a wide range of halo 
concentrations c = [5, 30] and disk mass fractions rn^ — 
[0.02,0.2]. However, such models often fail to reproduce 
the observed TF zero point (e.g. iNavarro fc Stcinmet j 



120001 iDutton et all 120071 ). Disk formation models 
can simultaneously fit the TF relation and the galaxy 
luminosity function if there was a ctually a sligh t halo 
expansion, where Kiax ^ 1-2 (jDutton et al.i r2007'). 
Evidence against adiabatic con traction is mounting also 
at z ~ 2 (jBurkert et al.ll2010t ). We find that a similar 
factor of about 1.2 between V^ax and Vh provides a good 
match to the TF zero point at z ~ 2. 

3.4. Redshift dependence 

According to Equation [TJ the observed SFR-M^ 
relation (equivalently the sSFR) at a given stellar mass) 
is declining with time from z = 2 to the present by a 
factor ^ 20. Not only does our model reproduce this 
evolution, but it can help us understand its origin. 

Figure [HJa) shows the model sSFR as a function 
of redshift for galaxies of different stellar masses, 
= 10^■^ 10^° and 10^°-^ Mq, in comparison 
with observe d sSFR for M^, i n the range 10^° — 
10^° -^ Mp, bylD addi ct al. (2009^); iDamen et afl ()2009aD 
and iPannellaet al. (2009 ) . One sees that our toy model 
reproduces the variation of sSFR with redshift relatively 
well to z = 2. We discuss the recent results at z = 4 to 
z — 7 in the Appendix. 

The redshift dependence of the sSFR is driven 
by redshift evolution of the average accretion rate. 
Equation \E\ which is oc (1 -|- z)^'^. Simply put, this 
factor is driven by the expansion of the universe. This 
can be clearly seen using the EPS formalism, where the 
factor (1 -I- z)^-^ originates directly from the derivative 
of the linear growth factor D ([Neistein fc Dekell [200l 
iBirnboim et al.l[2007| ) . Indeed, under the EPS formalism, 
the specific accretion rate is Mh/.^h = —s{M\i)uj, where 
s(Mh) is the self-invariant mass variable, and a;(z) = 
1.69/D{z) is the self-invariant time variable. s(Mh) is a 
function of the variance of the initial density fluctuations 
on mass scale My, . cr(Mh) (see Equation A5-A6 in 
IBirnboim et al.ll2007| ) 

The dotted line in Figure ^a.) demonstrates that 
the time evolution of the sSFR (Equation [6]) is indeed 
roughly proportional to uj, except that the predicted 
decline with time of the overall accretion rate is a bit 
slower than the decline in sSFR at z < 1 (Equation [1]) . 
An improved flt to the evolution at low redshifts is 
achieved by a small modification of the form uj x /(z) 
(long dashed line, defined in Section 12.31) . attempting 
to correct for the gradual decline with time of the 
cold gas fraction in the overall baryonic accretion. 
This confirms that a gradual decrease in the accretion 
efficiency parameter, by a factor of ~ 2 between z = 2 
and z = 0, helps recovering the observed evolution of the 
s SFR. 

iDavel ()2008[ ) parameterized the evolution of the SFR 
sequence at a given in terms of a star-formation 
activity parameter, defined as Ugi = {M^/Mi,)/{tH — 
IGyr). He pointed out that most models predict that 
galaxies have a growth rate comparable to tn , the Hubble 
time, i.e., with a constant ctsfr — 1, whereas observations 
reveal a growth of agf by a factor of ^ 3 from z = 2 to 
z = (see Fig[S|) . They used this apparent discrepancy to 
argue that the initial stellar mass function (IMF) must 
be evolving in time. Figure [B^b) shows the evolution 
of Usf predicted by our fiducial 'accFloor' model for 
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Fig. 6. — Left: Evolution of specific star formation rate (sSFR=SFR/Af^.) at a fixed from our fiducial model ('accFloor'). The 
dotted line shows the time dependence of the accretion rate through the function (ii oc (1 + z)^-^ (see text). The long-dashed line shows 
the slight modification of the introduced function /(z) (defined in Section l2.3l l. Right: Evolution of the star formation activity parameter 
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triangles) and lPannella et al.l II2009I ) (solid squares) for galaxies with log = 10-10.5. The shaded area shows the hydrosimulation results 
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Fig. 7. — Left: Evolution of the sSFR-M* relation. Right: Evolution of the main sequence. In both panels, the solid lines show the 
isochrones of the scaling relations at z = 4, 2.2, 1.2 and z = 0. The dotted lines show the individual tracks for the simulated halos. 



logM* — 10.0 and 10.7 Mq. An increase of agf by a 
factor of ^ 3 between z — 1 and z = 0, not far from the 
observational result, without any evolution in the IMF, is 
a general feature of the time dependence of the accretion 
efficiency. 

The SFR sequence and sSFR(M^) at any given redshift 
are simply isochrones in their respective parameter plane. 
Our model gives us the opportunity to show the actual 
tracks of individual galaxies, which are shown in Figure[7] 
for the sSFR-Af* (left) and the SFR-M^ (right) relations. 

4. ADDITIONAL ROLE OF M^,^ 



Given that our toy model seems to reproduce the 
scaling relations and their evolution rather well, we now 
compare the predictions of the model to observations 
involving the 'downsizing' effect (Section 14. ip . the halo 
stellar fractions (Section 14.21) . and the galaxy gas 
fractions fSection H31) . 

4.1. Downsizing and tau models 

Stellar-population analyses of red galaxies in the local 
universe show that stars tend to have formed earlier and 
over a shorter time-spa n in the m ost massive one s (e.g. 
iHeavens et al.l [2004 : .Thomas et al..,2005; Jimenez et al.l 
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Fig. 8. — Star- formation history SFR(t) for our simulated halos of 
decreasing mass (from left to right). The individual tracks simply 
reflect that more massive halos reach the mass floor M^[^ earlier 
as shown in Figure [T] and form stars on a shorter time scale (before 
reaching the mass ceiling Afmax). 



l2007HThoinas et al.ll2010[ ). This is sometimes termed as 
'archeological downsizing'. This age- mass correlation is 
also present in large samples (up to 3 x 10^) of active 
galaxies in the Sloan Digital Sky Surv ey (SDSS) (e.g. 
iGallazzi et al.ll2005l: [Jimenez et al.l l2007). 

The SF history (SFH) tracks shown in Figure M 
demonstrate the origin of the age-mass anti-correlation. 
More massive halos reach the mass floor earlier, and 
spend less time between the mass floor Mmin and the 
mass ceiling Mmax- This is a natural consequence of 
the halo tracks in Figure [T] combined with our key 
assumption of M^i^. 

Figure [9] shows the resulting age-mass relation at 
z — (left) and z — 2 (right). The left panel shows 
the formation redshift Zf as a function of stellar mass. 
Similarly to Noeske et al. (2007b), the formation redshift 
Zf is defined here when star formation began, i.e., when 
Mh{z) reaches Mmin- The mass index of the resulting 
age- mass relation is surprisingly close to the {1 + Zf) oc 
M°-3 inferred by INoeske et all ()2007b( ) shown as the 
dashed line (normalization arbitrary). 

That the SFR sequence leads to the same 
mass-dependent SFH as in the archeo l ogical d ownsizing 
was already argued by INoeske et al.l (|2007bD . Indeed, 
they fitted r-models to the z ~ 1 SFR-M^, relation and 
concluded the sSFR relation required both r and the 
formation redshift z/ to be mass-dependent (r oc M~^, 
(1 -I- Zf) oc M°-3). They cah this mass-dependent SFH 
'Staged Galaxy Formation' where more massive galaxies 
are formed first on a shorter time scale. However, the 
mass dependency in our model did no t have to agree 
with the resuhs of INoeske et al.l (|2007al ). 

This age-mass relation may already be in place at 
z ^ 2.2 (Figure Efrig ht)). Indeed, the large sa mple of 
the SINS survey (jForster Schreiber et all 120090 shows 
a similar trend between spectral energy distribution 
(SED) derived ages and stellar masses. In Figure [9l 
we compare these observations (gray points) to the 



luminosity weighted ages in our model (solid squares), 
which are derived according to / iibtSFR(tibt)dtibt where 
^ibt is the lookback time since z = 2.2. 

Thus, our model reveals that the key player for 
downsizing is the accretion mass floor Afmin. As it was 
first realized in .Neistein et al. (2006) , the archeological 
downsizing effect is a direct consequence of a mass fioor 

It has been noted that, for early type galaxies, the 
more massive ones reach the red sequence earl ier, a 
phen omenon dubbed 'downsizing in mass' (Cimatti et al.l 
l2006h . In our model, when galaxies reach the ceiling mass 
A^max, their gas supply dries out and SF is quenched, 
i.e., the galaxy becomes red and passive. The same 
SFH tracks shown in Figure |8] lead to this other type 
of downsizing, in t his case due to Mmax (see also 
iCattaneo et alJl2008| ). 



4.2. Stellar fractions 



shown that the stellar 
strong function of the 



fraction in 
total halo 



Many have 
galaxies i s a 

mass (e.e. Ivan d en Bo sch et al.l l2003bl: lEke et al.l 120051: 
iShankar et al."2006; van den Bosc h et al. l| 2007l : lKravtsovl 
2010; Mostcr ct al. 20101: iGuo eT al. 2010) owing to the 
very different shapes of the halo mass function and the 
galaxy stellar mass function. We confront our toy model 
with this additional observational constraint. 
Figure HOl compares the halo stellar fractions, defined 



as 



0.18 * Mh 



(15) 



predicted from our model to relevant data at z = 
(left) and z = 2.2 (right). The sha ded are a shows the 
results from the SDSS analysis of .Guo et ah! (l2010i). and 
the so lid line shows the similar result from lMoster et al.l 
(j2010t) . At z = 2.2, the observational estimates are 
from the disk-dominated galaxies of the SINS survey 
(jForster Schreiber et al.ll2009HCresci et al.ll2009D . where 
the virial mass is estimated assuming that = Vmax- 

We see that the models 'noMmin' and 'sfrFloor' 
overpredict the stellar fraction for halos below 10^^ Mq, 
while the 'accFloor' model recovers the shape /^.(Mh), 
including the drop toward lower masses occurring at 
the proper range of Afh- Once again, we learn that an 
effective halo mass fioor for accretion at Afmin ^ 10^^ Mq 
is essential for a fit to the data. 

However, the amplitude of the /*(Mh) function 
predicted in the simple accFloor model needs to be 
reduced roughly by a factor two for a match with the 
data. A better fit is obtained if we include in our model 
(Equation [3]) an outfiow term 

Afout -0.6M,. 

This constraint i s very con sistent with th e observations of 
IHeckman eFall (|2000l) and lMartinI H2005D who found that 

Aiout — SFR. Indeed, the corrections involved in deriving 
the outfiow rates in galaxies (for ionization, metallicity, 
and depletion) are usually of several orders of magnitude. 
Moreover, the instantaneous outfiow rate needs not be 
equal to our net outflow Afout term, which only describes 
the baryons that leave the halo forever. We note that 
this constraint on the outflow rate is degenerate with the 
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Fig. 9. — Archeological downsizing. Left: Formation redshift Zj (defined when M^^{z) reaches M^i^) as a function of stellar mass for 
our modeled galaxies (solid squares). The dashed line shows the Af* — zj relation inferred bv lNoeske et al.1 l l2007bl ) from the z ~ 1.0 SFR 
sequence, namely, (1 + Zf) oc M^'^ (normalization is arbitrary). Right: The modeled z ~ 2.2 age-Af* relation is shown with the solid 
squares. Modelled ages are the luminosity weighted ag es (see text). The gray point s show the ages (inferred from SED fitting) for the 
large sample of z = 2 SFGs available in the SINS survey l IForster Schreiber et al]|2009l ). At all epochs, the age- mass dependence is a direct 
consequence of the mass floor AI^i^. 
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Fig. 10. — Stellar fractions (/* = A/^/(0.18 * Afh)) at z = (left) and z = 2 (right) predicted from our fiducial model ('accFloor+'). In 
each panel, the alternative models, namely, with no mass floor Afnjin (dotted line), with an accre tion floor (long dashed line), and with an 
SF floor (short dashed line) are indicated. The 2 = stellar fractions from IMoster et al.l 12010) and IGuo et al.l 1120 10') are shown as the 
shaded area and thick line, respectively. The stellar fractions in z = 2 SFGs are shown as gray circles jf^rescT'et al.l l2009V 



(0.4—0.8) for ej, 



accretion efficiency ejn, i.e., 
ranging from 0.5 to 1.0. 

It is important to realize that the results presented 
thus far are independent of this outflow term. The SFR 
sequence is not sensitive to feedback as lowering the SFR 
by a factor of t wo leads to a s tellar mass Mi, lower by 
the same factor. iDutton et al.l (|2010( ) demonstrates this 
in a similar work. 

4.3. Gas fractions 



Another interesting constraint is provided by the gas 
fraction in galaxies, defined as 



/g 



+ Mi, 



(16) 



At l ow redshift, we learn from all-sky Hi surveys 
e.g. iZwaan fc et al.l l2003l: ' Rosenberg fc Schneideii 120031: 



Mever et al.l 



_ , , , . . _ _ g &: ocnneic 

20041: Giovan elU et all 120051 IWong eFall 
20061: iCatinella et al.ll20iai that the gas frac tions in SFGs 
vary systematically with stellar mass (e.g. iBaldrv et alj 
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120081 and references therein). These data are presented 
in the left panel of Figure [TTJ At z ^ 2, thanks 
to the rapid progress in the sensitivity of millimeter 
interferometers (such as the IRAM Plateau de Bure 
interferometer, PdBI), it is now possible to study CO 
rotational emission lines of normal high redshift SFGs 
at z > 1 and therefore constrain their gas masses. The 
observed gas fraction for z ~ 1-2. 5 SFGs with M^, > 
10^°-^M(7) ranges fr om 0.3 to 0.6 (jPaddi et all l2010al: 
ITacconi et al.l ^Olol ) with an a verage of about 0.44 for 
the large sample of 19 SFGs of ITacconi et~an (pOlOl ). 

Figure [TT] compares the predicted /g(M^,) from our 
model 'accFloor+' to those observations at z = (left) 
and z — 2 (right). We see that the model is consistent 
with the data both at z = and at z ~ 2. Locally, the 
upturn of fg at the low mass end is again a result of the 
mass floor Mmin- 

Contrary to the other results presented so far, the gas 
fractions are sensitive to SF efficiency ggfr, given that 
the SFR is set by the steady-state solution (or by the 
accretion rate) . A lower SF efficiency leads to higher gas 
fractions and vice-versa. We used a universal egfr set by 
the KS relation and the redshift evolution of fg in our 
model appears to be very consistent with observations. 
Indeed, the model predicts fg — 0.4 — 0.55 at z = 2 and 
fg — 0.3 — 0.45 at z — 1, which is in good agreement 
with the average gas fractions estimated from the CO 
obser vations, 0.44 and 0.34, respectively ()Tacconi et alj 
[20T0I see also Daddi et al. 2010). 

Our simplistic toy model, which was only tuned (with 
A'^min) to reproduce the scaling relations, helps us also 
to understand the trends of gas fraction with mass and 
redshift for SFGs. 



5. STAR-FORMATION HISTORY: CONSTRAINING 
THE MASS FLOOR 

Thus, based on the two scaling relations, the stellar and 
gas fractions, our model shows that a mass floor M^i^ 
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Fig. 12. — The star formation history p*. The dashed curves 
show that pir{z) is very sensitive to the minimum mass Mmj^. 
The light symbols show t he most recent z < 2 observations 
(Schim inovich et al.1 [20051 : IDahlen e t "all [20071: IMobasher et al 
2009) with the 2 > 2 data (Hopkins L Beacomll200g : IVerma et al 
2007; Wilkins et al. 2008) all converted to a Chabrier IMF 
following IHopkins fc BeacomI 1120081 '). 

must play an important role in galaxy formation and we 
argued that Mj^in should be '^10^^ Mq. We are about 
to show that the evolution of the average density of SF, 
/9^(z), will give us a quantitative constraint on Mmin- 

Figure [12] shows a compilation of data for the cosmic 
SFH, derived from observations under the assumption 
of a Chabrier IMF. The average density of SF, p*, is 
observed to be rising from z ~ 6 to z ~ 2, reaching a 
plateau at about z ^ 2, and then dropping by an order 
of magnitude between z ~ 1 and z = (e.g. iLillv et all 
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119961 iMadau et anil996t iHopkins fc Beacom|[200l and 
references therein). 

Our 'accFloor+' model convolved with the 
Sheth-Tormen mo del for the halo mass function 
(jSheth et alJ l2001[ ) is shown for 4 different values of 
Mmin, ranging from lO'^ to 3 x 10" Mq. We see that the 
value of Minin indeed determines the growth of p^, with 
time at high redshift. Best agreement to the growth rate 
observed at z > 2 is obtained with Mmin =10^^ Mq. 
Thus, we conclude that the rise in the cosmic SFH from 
z = 6 to z ~ 2 confirms the necessity of a mass fioor and 
provides the strongest constraint on the value of Mmin- 

Our model can also help us to understand the origin 
of the decline at late times. Possible reasons for this 
sharp decline could in principle be the acceleration of the 
universe at recent cosmological epochs, gas exhaustion 
into stars and outfiows, the quenching of SFR above 
the mass ceiling, or simply the generic decline in 
the accretion rate into dark-matter halos due to the 
expansion of the universe. The sharp decline at z < 2 is 
predominantly driven by the decreasing accretion rate as 
a direct result of the expansion of the universe as it can be 
seen from Figure [1] (right) and discussed in Section 13.41 
All the other potential sources for the decline, including 
the acceleration of the universe and the introduction of 
Mmax or gas consumption, have only minor effects on 
this decline. 

Ma ny other group s (e.g. iSpringel fc Hernquisti 
200l iNagamine et al.l [200a iSchave et al.l 120101: 
Choi fc Naeamind 120091) have used cosmological 
hydrosimulations to investigate the cosmic SFH 
and to identify the physical processes that drive its 
evolut ion. Similarly to our model, Hern quist fc Springell 
([2003) used an analytical approach and reach similar 
conclusions. They find that, at high-z, the rise in Pi,(t) 
originates from the gravitationally driven growth of 
halos, while at low-z, the decline is due to the expansion 
of the universe, through the inhibition of cooling. While 
in our model the SFR at a given halo mass is driven by 
the halo growth rate, their model assumes an ad-hoc 
form for the sSFR that is constant as a function of mass 
for all halos of virial temperature above 10"* K. Similarly 
to our A/min assumption, in their model, stars do not 
form in halos below Tvir = 10'* K. Their assumption also 
leads to a peak in the SFH, albeit at a higher redshift 
z ~ 5 or 6. In spite of the differences between their 
results and our results, this again illustrates the impact 
of a mass fioor. 

6. DISCUSSION AND CONCLUSIONS 

The purpose of this paper is to seek a basic 
understanding of the relations between global properties 
of SFGs and their time evolution between z ^ 2 and 
today. This includes the origin of the SFR sequence 
and the Tully-Fisher relation. We achieved this via 
an idealized and remarkably simple 'toy' model, based 
on the cosmological evolution of accretion rate into 
dark-matter halos, the efficient deep penetration of cold 
gas streams into the central galaxies and the standard 
efficiency of SFR in galactic disks. The model should be 
interpreted as a learning tool for gaining insight into the 
role played by several key physical processes. 

The model naturally leads to the observed scaling 



relations. Mi, oc M°-^ and Af* oc V^^x ^^'^ *o ^^i^ 
archeological 'downsizing' phenomenon, provided that 
there is a mass fioor Af,„in '^lO** Mq below which the 
accretion of cold gas is substantially suppressed 0. An 
analysis of the SFH allowed us to put an additional 
strong constraint on the mass fioor: Af,„in ~10^* Mq, 

corresponding to Vh ^ 100 A/hYf (1 + -^)3^2^km s^*, best 
accounts for t he observed sc al ing re lations. 

Recently, iDutton et al.l ()2010[ ) have used full 
semi-analytical approach to study the origin of the 
SFR sequence . Wh ile we agree with the conclusions of 
IDutton et al.l (|2010f ) with respect to the steady state, 
the impact of feedback, and to the evolution of the SFR 
sequence, our model differs from t heirs on one sign ificant 
aspect: the lack of a mass fioor. iCattaneo et al.l ()26l0f ) 
used a full semi-analytical model with an effective 
efficiency ein similar to ours (compare their Fig. 4 to 
our Fig.l). Their effective mass fioor is due to the 
combined effect of SN feedback and photo-heating 
due to reionization. These authors, however, used a 
re- ionization threshold at Vrc = 40 km s^* in order 
to fit the luminosity function, which is higher than 
expected (Kr ~ 10 km s~*) from pure r eionization (e.g. 
iBarkana fc Loebl 119991 : llhev et al.lT2005f ) by an order of 
magnitude in halo mass. 

6.1. Caveats and limitations 

Our model, which relies mainly on the cosmological 
accretion rate, is based on very few parameters (ein, 
Afmax, and Afmin) coui pared to semi-analytic models (e.g. 
ISomerville et"aLll20?)8h . Our main free parameter is A^min 
(affecting the cooling/accretion efficiency ein{M)), as we 
use normal a S F eff iciency e.sfr and set Afmax following 
ICattaneo et"all (e.g. l200l 12001 . 

The remarkably few paramet ers in our model are 
clearly an advantage (see also iNeistein fc WeinmannI 
l2009f ). but this inevitably leads to over-simplifications. 
For instance, modeling the lower and upper thresholds 
for cold gas accretion as sharp cutoffs is clearly 
unphysical. In reality, these transitions are likely to 
be smooth and to depend on other galaxy properties. 
For example. A/max is rather sensitive to the metallicity 
in the halo (Dekel fc Birnboim 200^). Indeed, the 
smoothness of this transition near A fmax can be derived 
from hydrodynamica l simulations ( Keres et al.l 120051: 
Birnboim etall l2007t lOcvirk et all boOSi: iKeres et al.l 
20091) . In fact, a smooth transition at Afmax would relieve 
the tension between the model and the observed scaling 
relations at the massive end, Af* ^ 10** Mq, in Figuresjl] 
and[5l 

The toy model presented in this paper does not allow 
us to address other observational constraints such as the 
mass-metallicity relation because this is rather sensitive 
to the exact treatment of the metal content of the stellar, 
SN outflows combined with the unknown fraction of ISM 
that is entrained in the onflow. 

We also cannot address the scatter in the scaling 
relations as the model deals only with the mean global 
properties. Variation s in the mass asse mbly history are 
to be expected (e.g. iGenel et al.l 120081) and will likely 

^ We note that the threshold could have been defined in terms 
of a constant virial velocity Vj, , and this would have lead to similar 
results. 
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contribute to the scatter in the scahng relations (e.g. 
iDutton et aT][2010l ). 

Another hmitation of our simphfied model is that we 
do not include bulge growth, which can be an additional 
drain for SF disks in our model. 

6.2. Conclusions 

In spite of its simplistic assumptions, our model 
captures many of the observed properties of SF galaxies 
and their evolution. This approach helps us gain insight 
into the following questions: 

• Why are the z — 2 SFRs so high? What sets the 
SFR sequence? Once the timescale for SF becomes 
shorter than the accretion timescale, galaxies settle in 
a quasi-steady state where the SFR is driven by the 
accretion rate (Figure [3]). As a consequence, the high 
SFRs for 2 = 2 SFGs are driven by the high accretion 
rate of cold baryons. The steady state also implies that 
the instantaneous SFR is comparable to the average SFR 
((SFR)) for a few Gyr from z ~ 4 to z ~ 1.5 (see 
Figure [ijb)). This is consistent with observatio ns from 
the SINS survey (jForster Schreiber et aD l2009f ) where 
the ratio (SFR) /SFR is found to be ~ 1.0 in z = 2 SFGs. 

• What is the origin of the scaling relations of SFGs ? 
We find that a single parameter, namely, the mass floor 
Afmin ^IO^^Mq for accretion, is key in generating the 
correct slopes of both the SFR sequence and the TF 
relation (Figs. H] and [5]). Thus, the two scaling relations 
are not independent relations, as we show in Section [3.31 

We emphasize that the SFR sequence and its evolution 
are insensitive to the values of the efficiency parameters 
ein and egfr, because changes in either of these efhciencies 
affect both M^, and Mi,, and shift galaxies along the SFR 
sequence. The same argument hol ds for the impact o f 
feedback on the SFR sequence (see IDutton et al.ll201ClD . 
This is the reason for the failure of any model that 
attempts to reproduce the tilts of the scaling relation by 
varying either the SF efficie ncy or the fe edback efficiency 
(see [Pavel 120081 : IDamen et al.l l2009aL and references 
therein) . 

• What sets the evolution of the SFR-M^, relation? 
The redshift evolution of the SFR sequence (or sSFR) 
is driven predominantly by the redshift evolution of the 
cosmological DM accretion rate. Both the SFR sequence 
and the sSFR at a fixed mass follow the cosmological 
decline of average accretion rate onto halos as a result 
of the expansion of the universe (Figure [S]). A slight 
improved match with the observations is obtained when 
we force the accretion efficiency to decrease with time at 
z < 2. 

• What is the origin of archeological downsizing? Why 
did more massive galaxies form their stars earlier 
and over a shorter duration (e.g. Heavens et al.l 120041: 
Thomas et "all 120051: iGallazzi et al]|2005t iJimenez et al.l 
2007[) ? An inspection of the halo trajectories in 
(Figure [T|) reveals that more massive halos reach the mass 
floor earlier, and spend less time between the mass floor 
and the mass ceiling (see also Figure [8]). It is the mass 
floor M,„in th at naturally pr o duces the age-mass relation 
(Figure [ni as lNeistein et al.l (l2006t ) already pointed out. 

•Is the observed low stellar fraction compatible with 
the high accretion efficiency of cold gas? Having a high 



accretion efficiencies (ein >50%) is entirely compatible 
with the observed stellar mass fractions that are 
low (< 50%), i.e., A ~ 0.20 (Figure ^ provided 
a net outflow rate Mout oc 0.6 Mi, is included. This 
is consistent with observations Mout — 1.0 M^, given 
that a significant fraction of the out-flowing gas may 
be returned into the interst ellar medium by z = (e.g. 
lOppenheimer & Davel l2008'): 

• Why are the gas fractions large at z > I? The 
observed high gas fractions at > 1 arise naturally from 
the high, mass-dependent accretion efficiency and the 
low star-formation efficiency. The SFR is set by the 
accretion rate in the steady-state solution, while the gas 
fraction is determined by the SFR efficiency Cgfr. The 
model predicts mean gas fractions of 30%-45% at z = 1.2 
and 40%-55% at z = 2.2 fFig ureITT]) in good agree ment 
with the recent obse rvations of iTacconi et al.l (|2010[ ) and 
iDaddi et all (|2010a[ ). 

• What is the origin of the shape of the SFH? The 
decline in the SFH from z = 2 to the present 
is primordially due to the decline in accretion rate 
associated with the cosmological expansion. The rise of 
Pi, from z = 6 to z = 2 necessitate a mass floor near 
M„,in - 10" M© (Figure ini). 

The successes of our model are evidence for the central 
role played by the key ingredient of our reservoir model, 
namely, the halo growth rate in the standard cosmology. 
While our model reveals the importance of a mass floor 
Mnun for accretion, it does not address its physical origin. 

7. ON THE ORIGIN OF THE MASS FLOOR 

The exact mechanism that suppresses or holds 
accretion below Mmin may be due either to photo-heating 
associated with the fresh UV background after 
reioiiization (e.g. iThoul fc Weinberg) 119961: 'Quin n et al. | 
1996tlBarkana fc Loeblll9 99: Gncdin 2000: Diikstra"e^L| 
200 l, or to SN feedb a ck (e . g. J ckel & Silk TTgseT 
[Dutton fc van den BoschI 120091 : lOp penheimer fc Davel 
l2008f ) 1^*1 Photo-heating associated with re-ionization 
increases the mean temperature of the IGM to ^ 10^ K. 
This reduces the cooling rate of the hot gas and prevents 
the assembly of new low-mass galaxies with Vh < 
10km s~^(or Mh ^ 10^ Mq), and those that have already 
collapsed a rc relatively quickly photo- evaporated (e.g. 
iBarkana fc L oeb 1999; Iliev et al.l l2005D . SN explosions, 
like photo-heating from rcionization, provide a negative 
feedback on SF. If present simultaneously, photo-heating 
and SN feedback effects can mutually a mplify each 
others ability to suppress the SFR f|Pieri fc Ma7tdll2007l : 
Pawlik_& Schayc 2009). Lastly, the formation of cold 
(r ~ 10-30 K) molecular gas (necessary for SF) requires 
specific physical conditions in ad dition to high gas 
densi ties, such as high gas pressure (jBlitz fc Rosolowskvl 
l200^ and pre-enrichment. Whether these conditions are 
met in low mass halos is an open question. Furthermore, 
radiative feedback effects from massive stars during 
intense SF episodes could play a crucial role in d isrupting 
the GMCs as proposed bv lMurrav et all (|2010[ ). 

Note, in the conte xt of momentum driven winds 
HOppenheimer He Dav^ I200SI 1. our mass floor corresponds to a 
threshold in the mass loading factor jj where ?j — ^ oo below Afmin- 
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It is possible that no single mechanism can produce 
such a mass floor. For instance, since simple arguments 
on SN feedback give A__oc M*/Mh cx Mh'^ with 
/3 = 2/3 (Dckcl & SillJ ll986( ). an additional physical 
process must be at play. The combination of SN 
feedback and a (very high) reionization threshold can 
produce a steeper f^, as shown i n Cattaneo et al. (2010). 
Simila rly, Ivan den Bosch et all (|2003ai) and iMo et al.l 
(j2005( ) showed that artificially induced re-heating of the 
IGM to T - 10^ K directly leads to a mass-dependent 
efficiency However, this effects alone also cannot 
account for th e steep mass dependence of (e.g. 
ILu fc Mol lpQOTi). 

Lastly, iCantalupol ()2010[ ) argued that the missing 
ingredient could be the effect of photoionization by 
local sources on the incoming cooling gas. He showed 
that soft X-ray generated by SF efficiently alters the 
ionization state of atoms, such as O. This effectively 
removes the main coolants and increases the cooling 
times by orders of magnitude, essentially stopping 
accretion, preferentially in low-mass halos. While this is 
an attractive mechanism, further work is needed in order 
to understand quantitatively whether photo-ionization 
from local sources are able to regulate the cosmological 
gas accretion rates. 



One way to gain insight into which of these possible 
mechanisms are at play is to constrain the transition 
in accretion efficiency near Mmin (ei„ oc MiJ'), which is 
related to the slope /? of the relation f^, oc M^/Mh oc Mh^ 
under the quasi-steady state solution with 7] ^ j3. Our 
modelling indicates that rj should be steeper than unity. 
Obser vationally, various groups have determined that (3 
is > 2 (Shanka'r et al."2006'; 'Baldrv et al."25ol|Kravtsov| 
[2010t fMoster et al...2010; .Guo et al...20ia) . 
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APPENDIX 

SPECIFIC SFR IN THE EARLY UNIVERSE 

Upon completion of this work, several groups have used deep Hubble Spac e Telescope data to their limi t and 
published constraints on the SFR sequence beyond z = 4 (|Stark et"ar 2009) , at z = 5 (lYabe et alJ l200^ , and 
z = 7 ()Gonzalez et al.ll2010[ ). As shown in Fig. [H some groups (.Stark e t all 120091: iGonzalez erahl I201Q) find 
that the SFR sequence do es no longer evolve, i.e., the sSFR re mains constant from z ~ 4 to z ~ 7 , whil e 
iSchaerer fc de B arro^ (I2010D find moderate and lYabe et all (I2009D strong evolut i on (s ee also iSawicki et al.l I2007D . 
The analysis of Schaere r fc de BarrosI ()2010[ ) includes the sample of IGonzalez et al.l (|2010D with the added WFC3 data 
from lOesch et al.l (j20ld|) . Given that it is a challenge to derive accurate (dust-corrected) SFR and stellar masses at 
those redshifts in a self-consistent way, it may be premature to view the (non-) evolution of iGonzalez et al.. (;2010l ) 
beyond z = 2 as definitely understood. 

Assumptions in the model 

Would the non-evolution of sSFR be real claimed by IGonzalez et al.l ()2010f ) , could the discre pant r esults of our 
model be due to some of our assumptions? Our model fails to indicate the turn over of Gonzale z et al.l (|2010[ ) in the 
evolution of sSFR(z) beyond z = 2 as shown in Figure [131 We turn to some of our assumptions in order to examine 
the robustness of this failed prediction. Could our model be modified in such a way to account for this apparent 
non-evolution of sSFR? Our main assumption, the accretion mass floor M„iin, only alters the mass index of the SFR 
sequence, but does not alter its redshift evolution. An assumption certai n to break dow n at z > 4 is our assumed 
constant recycling fraction R. Indeed, i? is a strong function of stellar ages (jBruzual fc Cha riot 2003) for ages < 1 Gyr. 
In order to understand how this affects our calculations in the first Gyr, we note that R affects the SF efficiency. The 
recycled fraction R is low at early times, meaning that the gas consumption is more efficient via the (1 — i?)SFR factor 
in EquationjH Since the SF efficiency is increased, the steady state solution is reached faster, i.e., sSFR(2) is going to 
be set by the redshift dependence of the accretion rate sooner. Therefore, a more realistic R{t) will not work in the 
right direction. 

At those early times, the main assumption to break down in our model is related to the steady-state solution. 
As detailed in Section l3Tl galaxies in our model are out of the steady-state solution beyond z > 5. Therefore, the 
assumption made for the SF timescale (or efficiency), namely that tsh scales with the halo dynamical time t^yn may 
break down: gas (or fuel for SF) is pouring at a rate faster than it can be consumed implying that the SF efficiency 
may be far from the local value. A more efficient gas consumption (shorter SF time scale) at early times means that 
the steady-state solution is reached more rapidly. This will only exacerbate the discrepancy between the z =4-8 data 
and the model. A less efficient gas consumption (longer SF time scale) at early times means that the SFR is increasing 
to some power of time t or redshift (1 -I- z). Thus, the sSFR is always going to be proportional to sSFR(t) oc 1/t. 

Therefore, upon examining our assumptions, we find that the non-evolution of sSFR beyond z = 2 cannot be achieved 
easily. In other words, the evolution o f sSFR beyond z = 4 is giving strong constraints on early galaxy formation. 
Taken at face- value, the IGonzalez et all ()2010[ ) results are going to be extremely challenging for any theoretical model. 

Assumptions in the observations 

In the observati ons, in order to fit th e observed SEDs, one has to assume (i) a SFH, and (ii) treat the nebular emission 
lines consistently (jSchaerer fc de B"aiT OS 2009) . These two can affect the sSFR significantly. For instance, at high SFRs 
and young ages, the nebular emission lines (e.g., [O III], Hf3) can affect the color significantly (jSchaerer fc de BarrosI 
11009; Yabe et al. 2009, Finlator, in prep.). When this is omitted, the resulting Balmer break is overestimated, and 
hence the stellar mass, leading to an underestimated sSFR. 

Similarly, the SF H can have a strong impact on the derived sSFR. At z = 2, the often assumed 'constant SFH' works 
generally well (e.g. iForster Schreiber et al.ll2009l) and this assumption is in good agreement wit h the growth predicted 
(see Fig.|6|and Section IO]) . At z > 5, diff e rent g roups make different assumptions. For example. IGonz alez et al.l ()2010f ) 
prefer a constant SFH, while IStark et al.l (j2009( ) and Schaerer fc de Barros (2010) use an exponentially declining SFH 
with a range of ages. Unfortunately, the SFH is difficult to constrain from the data themselves, but has a significant 
impact on the sSFR. Indeed, if SFR(f;) is constant, then sSFR goes as 1/t. If SFR(t) is declining exponentially, then 
Mi,{t) is dominated by the initial SFR at t = to {Mi, ~ exp(— to/T")), and sSFR is dominated by the age At = t — to, 
with sSFRoc exp(— At/r). If SFR(t) is increasing exponentially, then M^(t) is proportional to exp(-|-t/T), and the 
current SFR is also cx exp(-|-t/r). Hence, the sSFR is constant in this case. 

Given the short times between z = 5, z = 6, and z = 7, the various populations observed at these epochs are likely 
descendants and progenitors of one another and the inferred involution of sSFR(z) ought to be consistent with the 
assumed SFH. Currently, it appears that this self-consistency in the observational results is not present. 

Star formation histories 

At z > 5, galaxy formation is surely more uncertain. Proto-galaxies are more metal poor, the universe was just 
re-ionized, two facts that will impact molecular gas formation, hence the KS relation itself. Hence, several of our 
assumptions are likely to fail at those early epochs. However, a robust and generic feature of any model at those 
epochs is the rapid growth rate Af , which implies that the SFHs are very far from exponentially declining. In order to 
guide interested readers with SFHs appropriate in this regime of steady growth (at z > 4), we note that in the redshift 
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Fig. 13. — The evolution of sSFR at high-redshifts. The lines show the predicted predicted evolution by our model ('accFloor+') at a 
fixed Af* as in Fig. (6] At z < 4, the da ta points are as in Fig. E] At 2 > 4, the Gonzalez ct al. (2010) and Stark et al. (2009) results 
indicate no-evolution of the sSFR, while lYabe et al.l (12009 1) indicates strong evolution. This clear lack of agreement among the various 
groups refiects the challenge s in deriving accurate (dust-corrected) SFR and stellar masses at those redshifts in a self-consistent way (see 
text). The re-analysis of the I Gonzalez et al.l (|2010l 'l sample bv iSchaerer fcT le Barro^ (|2010| ) using the rece nt WFC3 data is shown with the 
open square with solid error bars (with logM* ~ 9.5). The 'faint' sample of l^chaerer &: de Barroj 120101) with logM* ~ 8.0 is shown with 
the square with dotted error bars. 

range z = 5 to z = 7, our SFHs are approximate increasing exponential, SFR(t)cx exp{+t / t) , with a time scale r close 
to ~ 0.5 Gyr as shown in Fig. [T3] with the dotted line. Using r = 0.5 Gyr, one can approximate the SFR(t) with 



SFR(i) cxe(+*/^)for times t > tf, 



(Al) 



where the formation time tj is given by (1 + Zf) ~ 0.4 x M'^-^^. As Fig. [Til shows that r varies from 0.3 Gyr to 0.6 Gyr 
between z = 5 and z = 7, a physically motivated form (dashed line) is 



SFR(t) = 100e(+(*-*^)/^(^»with t{z) = 0.5Gyr(it^'l 



+ Zb J 



(A2) 



where tc and Zf, are the fitted parameters. Table [2] lists the fitted parameters. A more accurate fit is obtained using 
polynomials of the type gq + aiU^ + 02?/^, where U = t/ (lO^yr), and shown as the red curves in Fig. [TH 
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Fig. 14. — At z > 4.5, the star-formation histories are close to being increasing exponential, SFR(t) oc e"'"*/'^ with r ~ 0.5 Gyr (dotted 

/ . \ -2.2 

line). Motivated by the halo growth, the dashed line shows an evolving tau r{z) ~ 0.5Gyr( ^"^^^ j . The red curves shows second order 
polynomial fits. 
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TABLE 2 

Parameters for the SFHs using Equation IA2I . 



